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Abstract 
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1  Introduction 

The  linear  time  varying  descriptor  system 

=  m  a) 

plays  an  important  role  in  the  growing  understanding  of  descriptor  systems. 
It  is  a  bridge  between  the  well  understood  linear  time  invariant  problem 

Ex'(t)  +  F  x(t)  =  f(t)  (2) 

and  the  general  nonlinear  implicit  differential  equation 

G(x',x,t)  =  0  (3) 

Many  of  the  difficulties  that  separate  (2)  from  (3)  already  occur  in  the  linear 
time  varying  case  (1).  On  the  other  hand,  the  linearity  of  (1)  makes  possible 
a  more  complete  analysis.  Understanding  developed  from  analyzing  (1)  has 
proved  useful  in  developing  numerical  and  analytic  techniques  for  (3)  [6]. 

The  first  general  numerical  method  for  (1)  is  developed  in  [2], [3].  This 
method  is  based  on  working  with  arrays  of  derivatives  of  E,  F,  f.  It  has  the 
advantage  that  all  differentiations  are  performed  directly  on  the  coefficients 
of  (1)  without  any  time  varying  coordinate  changes.  Theoretical  character¬ 
izations  of  solvability  [4]  for  (1)  and  a  numerical  procedure  for  computing 
a  state  space  realization  [5]  have  come  out  of  this  approach. 

Simultaneously,  there  has  been  progress  in  studying  implicit  differential 
equations  using  the  ideas  of  geometric  control  theory.  In  particular,  [11] 
considers  realizations  utilizing  Lie  derivatives. 

Intuitively,  there  is  a  close  connection  between  realizations  and  numeri¬ 
cal  methods.  Many  numerical  ordinary  differential  equation  (ODE)  integra¬ 
tors  are  based  on  being  able  to  estimate  x'  given  x,t,  whereas  a  realization 
is  often  in  the  form  of  an  ordinary  differential  equation  x'  =  Q{x,  t). 

This  paper  is  the  beginning  of  our  effort  to  unify  these  two  approaches. 
There  are  several  potential  benefits  to  this  effort.  Hopefully,  this  will  make 
the  tools  of  geometric  control  theory  available  in  trying  to  develop  and  ana¬ 
lyze  numerical  methods  for  •:!)  t'cnversely.  the  results  from  the  numerical 
theory  will  suggest  better  a> «  to  compute  the  geometric  objects.  Also, 
as  will  be  shown  in  thia  pap-r  c;i:ic  ideas  from  the  numerical  theory,  we 
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will  be  able  to  extend  some  of  the  geometric  theory  thereby  obtaining  new 
theoretical  results. 

Section  2  will  quickly  summarize  the  geometric  method  of  realization 
given  in  [11].  Section  3  summarizes  the  key  ideas  from  [4].  In  Section 
4,  we  develop  the  relationship  between  the  ideas  in  Sections  2  and  3  and 
establish  extensions  to  the  theory.  Finally,  Section  5  contains  concluding 
remarks  and  discussion. 

2  Geometric  Realizations 

This  section  will  give  an  outline  of  part  of  the  method  of  realizations  for  (1) 
developed  in  [11].  An  exposition  of  the  underlying  nonlinear  control  theory 
can  be  found  in  [9]. 

Consider  a  nonlinear  state  system,  x  ETln  \ 

x'  =0o(*)  +  ]E>J^(*)  (4) 

;=i 

A  submanifold  M  C  TV  is  controlled  invariant  if  there  exists  a  smooth 
feedback  =  a7(x)  such  that  when  applied  to  (4)  it  makes  the  vector  field 
for  (4)  tangent  to  N.  That  is,  solutions  starting  in  A/",  stay  in  A/*.  For 
a  scalar  function  T(x)  =  T(x i, . . .  ,  x„)  from  TV  to  Tl,  and  vector  valued 
function  s(x)  from  TV  to  TV,  let  L,T  be  the  derivative  of  T  along  the 
differential  equation  x'  =  s(x).  Thus 

d  "  dT  dx,  ^  dT  ,  .  __  .  , 

L.T  =jtT  =  T. tetlr  =  E £-*<*>  =  vr  .(*) 

Notice  that  L,T  is  again  a  scalar  valued  function  of  x  so  that  L]T  = 
L,(L,T)  is  well  defined  and  similarly  L\T  is  well  defined  for  nonnegative 
integers  r. 

We  are  interested  in  the  case  when  (4)  satisfies  constraints 

H(x)  =  0  (5) 

where  H  has  values  in  TZP.  As  in  [11]  consider  p  output  equations 

y  =  H{x)  (6) 
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Let  yi  =  Hi(x)  be  the  ith  equation  in  (6).  For  i  =  1  ,  ...,p,  let  p,  be 
the  smallest  integer  p  such  that  L^L^/f^x)  ^  0  for  some  x  and  some 
j  €  {1, . . . ,  m}  Assume,  as  in  [11],  that  p,  exists  for  i  =  1, . . .  ,p.  Now  let 
the  p  x  m  matrix  A(x)  be  defined  by 


AtJ  =  L3j  Lg'0H,(x) 


(<) 


The  next  assumption  in  [11]  is  that  A  has  constant  rank  p.  Then  the 
maximal  controlled  invariant  submanifold  of  H(x)  =  0  exists  and  is  given 
by 


AT  =  {x  |  L^,(x)  =  0,  fc  =  0,...,p„  i  p} 


(8) 


Furthermore  the  needed  feedback  is  computed  as  follows.  Let  a(x)  satisfy 


A(x)q(x)  +  b(x)  =  7(x) 


(9) 


where 


b,(x)  =  L*+' H,(x)  (10) 

and  7  is  a  vector  of  functions  which  are  zero  on  jVm .  We  shall  take  7  =  0. 
Let  0(x)  bemx(m-p)  with  rank()3)  =  m - p such  that  A3  =  0.  Then  the 
feedback  law  u  =  a  +  3v  applied  to  (4)  gives  a  control  system  with  m  —  p 
inputs  v  which  leaves  Xm  invariant  and  thus  can  be  reduced  to  a  control 
system  on  Xm : 


m  — p 


X  =  go{x)  +  52  as(x)9j{x)  +  12  12 

j  =  l  k=l  \;=l  ) 


(ID 


Notice  that  this  theory  provides  two  key  objects.  The  manifold  A’* 
consists  of  solutions  of  (4)  that  satisfy  the  constraints  H[x)  =  0  for  some 
control  u  while  (11)  gives  a  control  system,  defined  on  the  whole  space, 
which  leaves  A/**  invariant  and  on  Xm  agrees  with  the  solutions  of  (4). (5). 

This  theory  is  applied  to  implicit  differential  equations  in  [11]  as  follows. 
Suppose  that  we  have  the  kth  order  descriptor  system. 


Rt[u>,w . u,ifc|]  =  0.  w  €  TV ,  i  =!,...,£ 


(12) 


I 

A. « 


i 

i'l 


...  /v 


This  is  rewritten  as 


j_  u;(fc> 


with  output 


0/0 


•  *  0 

•  0  I 

■  •  0 


z  =s  R[w,  w' , . .  . , 


Here  the  derivatives  of  w  are  considered  to  be  variables.  Note  that  w 
is  a  solution  of  (12)  if  and  only  if  there  is  a  v  so  that  the  solution  of 
(13)  makes  the  output  (14)  zero.  Writing  (13), (14)  as  (4), (6)  using  x  = 
(w,  tv', . . . ,  ui(fc)),  we  have 


g0(z)  = 


.  9j(*)  = 


where  {e^}  is  the  standard  basis  for  7 ln ,  and 


Then 


H(x)  =  /2[u»,  w  , . . . ,  u;(fc)] 


Aix)  = 


Assuming  that  the  rank  of  A(x)  is  /,  there  then  exists  a  feedback  u  = 
a(x)  +  0(x)v  so  that  x'  =  g0(x)  +  E  u;?;(x)  leaves  N'm  controlled  invariant 
and  hence  will  generate  a  realization  in  yV'*  coordinates. 

We  shall  illustrate  these  ideas  with  an  example  of  a  simple  singular 
system.  This  example  will  be  referred  to  again  later. 

Example  1  Consider  (2)  with 


-1  0 
0  -1 


j-1  itwjtkjtw  u 


so  that 


Zj  =  Zi  +  u, 
0  =  X'l  +  U2 


and  /  =  u.  Letting  vj  —  [x,u]  =  [xt,  x2,  Uj,  u2],  w'  =  [x',u' ]  we  get  that 
(13), (14)  are 

d[u/l  [0  J  1  [  u;  1  f  0  1  /10, 

dt  w'  ~  [  0  0  j  w'  +  [  /  U  (18^ 

and 


R\(w,  w')  =  X2  -  Z;  -  U\ 
R,2(w,w')  =  X2  +  U2 


t  w'  0  . 

9o{w,w)  =  0  ,  g}  =  ^  ,  ;  =  1,...,4 


To  simplify  our  notation,  let  0r  be  an  r  dimensional  zero  vector  and  let  ; 
separate  the  blocks  of  n  entries.  Then 


L^R,  =  L3lR\  =  [-1  0  -  1  0;0  1  0  0][04;ej]r  =  62} 
Hence  p\  —  0  since  L3iR\  =  1  ^  0.  However, 

L3>R2  =  [0  10  l;04][04;ejf  =  0 


£<70-^2  =  [010  1:  04][z'1 ,  x2,  u'j,  i4;  04]r  =  x'2  +  u2 


Hence 


I3,L3o/?2  =  [04:0  10  1][0;eJ']T 
is  1  if  j  =  2,4  and  zero  otherwise  Hence  P2  —  1  and 


[0  1  0  O' 
,0101 


Furthermore  X’  is  giv-n  (•> 


=  x\  -  Xi  —  Ui  =  0 
=  X;  —  U2  =  0 
=  x.  —  <u  =  0 


*  Va'.V.s’a  lV.v.V/. 
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or  equivalently, 


-«i  “  A 
-u2 


XX  - 
X]  = 


(22) 

(23) 


which  is  exactly  the  solution  manifold  of  (2)  for  this  example.  Continuing, 
we  compute  that 


Let  7  =  0.  Then  one  solution  for  a,/?  is 


0 

'  1 

0  ' 

x\  +  u\ 

II 

0 

0 

0 

0 

1 

l - 

l 

-N- 

1 

JP- 

i _ 

0 

0 

Applying  the  feedback  v  =  a  +  gives  the  explicit  control  system  (11) 
which  after  deleting  the  redundant  equations  is 


x" 

= 

t>i 

x"i 

— 

Xj  -1-  u 

n 

_ 

U1 

= 

t’2 

ii 

“2 

— 

“*1  ~ 

(24) 


This  explicit  system,  for  the  correct  choices  of  vx,  v-2  includes  the  solutions 
(22), (23)  of  the  original  implicit  descriptor  system  if  the  initial  conditions 
satisfy  the  jVm  conditions  (21). 


3  Derivative  Arrays 

As  in  [4]  we  assume  that  E .  F  f  are  infinitely  differentiable  to  avoid  tech¬ 
nical  difficulties.  The  system  i  1 1  is  solvable  on  an  interval  1  if  for  every 
/,  there  is  a  solution  x  < on  all  of  1  and  for  a  given  /,  solutions  are 
uniquely  determined  by  'h~:r  value  at  any  t0  €  1.  See  [4]  for  a  more  careful 


definition.  Let  for  c  —  E,F,f,x.  Then  for  any  j  >  0; 

£^°1  0  •  0  1  *M  ' 

£W  +  ir(0]  2E[a\ 

*  *  *  o 

£-h-i)  +  fb- 2]  2 +  i*~3S  *  j E[0]  *M 


'  /( 0]  ■ 

jrfo] 

fb- 1] 

/rb-i) 

where  =  x.  We  rewrite  (25)  as 

£jXj  =  fj-FjiW  (26) 

The  matrix  Ej  is  said  to  be  1-full  if  there  is  a  nonsingular  matrix  ©  such 
that 

_  /  0 

e€]  ~  0  M 

where  I  is  n  x  n.  From  [4]  we  have  the  following  fundamental  result; 

Theorem  3.1  The  system  (1)  with  E,F  infinitely  differentiable,  is  solvable 
if  and  only  if  there  is  a  j  such  that 

(i)  Ej  has  constant  rank  on  T 

(ii)  Ej  is  1-full  for  every  t  e  I 

(ui)  7 =  Cn  for  all  t^l,  1  <i<j 

Furthermore,  tf  (i)-(iii)  hold  for  j0,  then  they  hold  for  any  j  >  j0.  Also, 
j  —  n  +  1  satisfies  (i)-(iu). 

The  result  we  shall  need  from  [3], [4]  is; 

Theorem  3.2  Suppose  that  (1)  is  solvable  and  j  satisfies  (i)-(iu)  of  The¬ 
orem  1.  Then  the  manifold  of  consistent  initial  conditions  at  time  t  is 
characterized  by 

Ej(t)Xj  =  f fit)  -  Fj(t)x(t)  (27) 

is  consistent  That  is,  f,  —  F,x  el Z(E}). 


* 


A  numerical  method  utilizing  these  Theorems  is  developed  in  [3].  But 
for  our  purposes  here,  it  suffices  to  recall  that  we  can  multiply  (26)  by  a 
nonsingular  @  to  get 


'In  0  O' 

"  9i(0  ‘ 

'  Qi(t) ' 

0  MM 

X;  = 

<k(i) 

— 

Qi{i) 

.00  0 

.  (0  . 

[  (hit)  J 

The  first  row  of  (28)  is 

x'  =  -<2i(*)*  +  ?i(<)  (29) 

which  is  an  explicit  ordinary  differential  equation  whose  solutions  include 
the  solutions  of  the  original  implicit  differential  equation  (1).  The  third 
row  of  (28), 

Qs(t)x  =  q3(t)  .  (30) 

characterizes  the  solution  manifold. 


4  Discussion 


We  begin  our  discussion  by  considering  the  simple  solvable  system  in  Ex¬ 
ample  2. 

Example  2.  In  (2),  let 


1  - 
1  -1  I  ’ 


0 

1 


If  the  procedure  in  Section  2  is  applied,  we  get  p\  =  p2  =  0  and 


1-10  0 
1-10  0 


which  does  not  have  full  row  rank. 

Example  2  shows  that  the  procedure  of  Section  2  can  experience  diffi¬ 
culty  with  solvable  systems  Of  course,  for  this  example,  a  constant  coor¬ 
dinate  change  would  convert  it  to  Example  1.  However,  the  general  case 


0 


is  much  more  complicated.  There  may  not  exist  any  smooth  coordinate 
changes  that  put  a  solvable  system  into  a  nice  structural  form.  Even  if 
there  are  such  coordinate  changes,  they  will  be  time  varying,  even  nonlin¬ 
ear  in  the  case  of  (4),  and  thus  sometimes  difficult  to  compute.  In  general 
then,  if  time  varying  coordinate  changes  are  not  computed  and  the  method 
of  Section  2  is  being  used,  we  would  expect  the  matrix  A  to  have  varying 
rank  or  be  rank  deficient,  unless  the  system  (1)  has  a  well  defined  structural 
index. [8]  This  does  happen  with  some  systems  of  physical  importance,  such 
as  some  of  those  in  mechanics  [1],[7],[10]  but  is  not  a  general  characteristic 
of  solvable  systems. 

The  approach  of  Section  3,  like  that  of  Section  2,  produces  a  description 
of  the  solution  manifold  (30),  and  a  differential  equation  defined  on  all  of 
7Zn  (29)  whose  solutions  include  those  of  the  singular  system.  However,  in 
Section  3,  no  attempt  is  made  to  take  the  minimum  number  of  differenti¬ 
ations.  Rather  a  sufficient  number  are  taken.  Then  all  needed  algebra  is 
performed  pointwise  so  that  it  is  never  necessary  to  differentiate  computed 
quantities  which  is  numerically  unstable. 

In  this  section,  we  shall  discuss  the  adaptation  of  these  ideas  to  the 
approach  of  Section  2.  Our  long  range  goal  is  two  fold.  First,  we  want 
to  extend  the  theory  of  Section  2.  Secondly,  we  wish  to  utilize  the  insight 
gained  from  the  method  of  Section  2,  in  order  to  extend  the  approach 
of  Section  3  to  general  nonlinear  descriptor  systems.  Any  approach  which 
generates  an  explicit  differential  equation  from  an  implicit  one  is  potentially 
useful  in  making  possible  the  application  of  explicit  differential  equation 
numerical  integration  schemes  to  implicit  differential  equations. 

Consider  again  the  linear  time  varying  singular  system  (1).  Differenti¬ 
ating  the  equation  j  —  1  times  yields 


E 

E'  +  F 
E"  +  2  F' 

* 

+  (j  -  l)F,J-2) 


0 

E 

2  E'  +  F  E 


where  the  ij-entry  of  (31)  for  i  >  j  is  ^  _  “  J  EK'  ^  '  .  ~  J  F[l  J~l> . 

We  write  (31)  as 

ijXj  =  -?jX  +  fj  (32) 

Let  Dj  =Diag{l,  1,  ,  ^rijr},  D:  =Diag{l!,  2!, . . .  j!}.  Then 

DjTj  —  Fj,  Djij  =  f j ,  DjXj  =  Xj,  D}SjDj  —  S: 

Thus  Theorems  3. 1,3. 2  hold  using  (32)  instead  of  (26). 

Now  consider  (1).  Again  let  u  —  f  and  w  =  [z,u],  and  rewrite  (1)  as 

r'  =  1  (33) 

R(r,w,w')  —  E(r)x'  +  F(r)x  —  u  =  0  (34) 

If  we  were  following  Section  2  we  would  let  k  =  1  since  there  are  only  first 
derivatives  in  the  R  equation.  However,  we  know  that  in  general  additional 
derivatives  are  needed.  Consider  then  the  system 


W  2(j+l)n  r  n 

.  V'  Ul 

...in  «=i  e' 


2  =  R[r.u.  .  .  =  E{r)x'  +  F{t)x  -  u  =  0  (36) 

'  1  1 

'i  '  I  r 


'  xW.VJVlV 


Then 


is  the  same  as 


CD 

II 

cs 

Lg0  R  —  0 

(37) 

tr-< 

°  l 

II 

O 

(i)  +  Fj(r)x(t)  -  Uj  =  0 

(38) 

Thus  Theorems  3.1,  3.2  show  that  the  extended  family  of  Lie  derivatives 
(37)  will  determine  asolution  manifold  and  aflow  in  all  of7£n.  However,  this 
fact  is  expressed  in  terms  of  Section  3.  We  wish  to  exploit  the  information 
in  these  theorems  to  extend  the  approach  of  Section  2. 

If  a  solution  of  (35)  is  to  leave  (36)  invariant  for  a  control  v,  we  need 
that 

jt(LTgoR )  =  L'+'R  +  (LsLrgoR)v  =  0  (39) 

for  r  =  0, . . . ,  j  —  1,  (j  <  n)  where  the  ith  column  of  LgLrgoR  is  LgtLrgoR. 
The  equations  (39)  may  be  rewritten  sis 


b  +  Aa  =  0 


(40) 


and  A  =  \Tg ,  S3 ,  —  I]P  where  P  is  the  permutation  matrix  that  lists  the  x(,) 
first  and  then  the  u(,).  Notice  that  A  has  full  row  rank  by  Theorem  3.1.  As 
before,  define  the  feedback  i;  =  d  +  3v,  by  Aa  =  —b,Afi  =  0.  Using  these 
theorems,  observe  that  not  only  is  there  such  a  feedback,  but  that  in  3  we 
may  take  u  arbitrary  and  then  x  is  determined  by  u;,x0- 
Alternatively,  we  may  let  u-  =  x,  and  rewrite  (1)  as 


d_ 

dt 


x 


z  —  i?[r.  x. 


L  0  J 


(;  +  l)n 

£  ‘v 


i=i 


nl  r)x'  —  F{r)x  —  u(r)  =  0 


(41) 


(42) 
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Again  (37)  becomes  (38).  This  time  A  =  [Tj,Sf\  and  the  equations  (39) 
are 

b  +  it,  =  [Tj,£j]6i  (43) 

From  Theorem  3.1  we  again  have  that  A  has  full  row  rank  and  hence 
d  exists.  Furthermore,  in  this  case,  the  second  n  components  of  /?  depend 
only  on  the  first  n  components  of  /3 

Another  variation  of  the  preceding  approach  would  be  to  consider 

E(t)x'  4-  F(t)x  =  B(t)u  (44) 

In  (35), (36)  and  (41), (42),  we  have  introduced  the  variable  r  to  make 
the  resulting  systems  time  invariant,  and  thus  more  closely  resemble  those 
in  Section  2.  However,  this  is  not  necessary. 

Which  of  these  alternatives  is  the  best  way  to  view  (1)  remains  to  be 
determined. 

5  Conclusion 

In  this  paper  we  have  briefly  discussed  two  approaches  concerning  descrip¬ 
tor  systems  and  observed  that  they  deal  with  similar  derivative  arrays. 
The  approach  of  Section  2  is  part  of  an  evolving  elegant  differential  geo¬ 
metric  control  theory.  The  approach  of  Section  3  is  computationally  simpler 
in  that  the  only  needed  symbolic  manipulations  are  derivatives  of  known 
functions.  Numerical  linear  algebra  routines  may  then  be  used  pointwise. 
The  assumptions  of  Section  3  are  also  somewhat  more  general  than  those 
of  Section  2. 

In  subsequent  papers,  we  hope  to  exploit  this  relationship  more  fully, 
especially  in  the  nonlinear  case. 
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